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The  Shape  of  a  Liquid  Drop  in  the  Flow  of  a  Perfect  Fluid 
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A  method  of  analysis  previously  used  to  determine  the  shape  of  a  liquid 
drop  in  an  electric  field  has  been  applied  to  determine  the  shape  of  a  drop  in 
the  flow  of  a  perfect  fluid.  The  surface  of  the  drop  is  expanded  in  a  series 
of  Legendre  polynomials  with  arbitrary  coefficients.  Minimisation  of  the 
total  energy  with  respect  to  the  coefficients  then  gives  equations  for  their 
determination.  Due  to  the  similarity  of  the  equations  for  determining  the 
coefficients  to  previous  results,  it  is  expected  that  the-  shape  of  the  drops 
will  deviate  considerably  from  the  usually  assumed  oblate  spheroids. 
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1 •  INTRODUCTION 

In  this  report  we  determine  the  shape  of  a  liquid  drop  in  the  flow  of  an 
ideal  fluid  (no  viscosity).  The  drop  is  assumed  to  be  held  together  by  the 
interface  surface  tension  between  the  liquid  inside  the  drop  and  the  external 
fluid.  In  previous  analyses  of  this  problem,  the  shapes  of  the  drops  have 
been  assumed  to  be  oblate  spheroids  with  a  minor  axis  parallel  to  the 
direction  of  flow  of  the  external  fluid. In  our  earlier  analysis  of  the 
shape  of  liquid  drops  in  an  electric  field,  the  commonly  accepted  assumption 
that  the  drops  formed  prolate  spheroids  with  the  major  axis  parallel  to  the 
electric  field  was  shown  to  be  generally  invalid.3'4  The  spheroidal 
approximation  for  drops  was  only  vail  *  for  small  drops  or  law  electric  fields. 

Since  the  spheroidal  approximation  was  invalid  for  drops  in  an  electric 
field,  we  suspected  that  perhaps  the  same  approximation  would  be  invalid  for 
falling  liquid  drops  in  air  or  other  fluids.  We  have  made  the  assumption  that 
the  external  fluid  has  no  viscosity;  otherwise,  the  problem  is  extremely 
complicated.  Here  we  shall  use  the  technique  used  in  our  previous 
analysis.4  That  is,  we  express  the  total  energy  of  the  system,  the  kinetic 
energy  of  the  external  fluid,  and  the  potential  energy  due  to  surface  tension, 
in  terms  of  the  coefficients,  aR,  in  the  expansion 


9 


r(8)  «  a  +  T  a  P  (cos  8) 
o  **  n  n 


(i) 


which  relates  an  arbitrary  point  on  the  surface  of  the  drop  to  the  angle  8 

measured  from  the  direction  of  the  fluid  flaw  (z  axis).  In  equation  (1)  the 

P  (cos  8)  are  Legendre  polynomials.  Minimizing  the  total  energy  with  respect 
n 

to  the  an's  gives  a  set  of  equations  for  these  coefficients.  The  solutions  of 
this  set  of  equations  then  are  used  in  equation  (1)  to  determine  the  shape  of 
the  drop. 


2.  THEORY 

In  equation  (1),  me  have  made  the  physically  realistic  assumption  that 
the  drop  is  symmetric  about  the  direction  of  the  external  fluid  flow.  In 
standard  texts  on  fluid  mechanics. 5,6  it  is  shown  that  the  velocity  in  a 
perfect  incompressible  fluid  can  be  derived  from  a  velocity 
potential.  +.  which  is  a  solution  of  Laplace's  equation.  That  is. 


and 


V2*  -  0  (2) 

▼  -  -V*  .  (3) 


The  solution  to  equation  (2)  appropriate  for  the  region  external  to  the  drop 
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where  we  have  dropped  the  arguments  of  the  Legendre  polynomials 

[p  ■  P  (cos  9)1  and  assumed  that  far  from  the  drop  the  velocity  of  the 
n  n 

external  fluid  is  given  by 


▼  -  eEV  ,  (5) 

A 

and  ez  is  a  unit  vector  in  the  s  direction. 

As  in  our  previous  analyses3'4  me  shall  rewrite  equation  (1)  as 

r(9)  -  a  ♦  8  7  a  P  ,  (6) 

o  **  n  n 


where  4  is  an  order  parameter.  We  then  expand  the  An  in  equation  (4)  in 
powers  of  4  as 
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(7) 


2 

and  retain  the  terms  through  4  . 

The  A_  of  equation  (4)  are  determined  by  the  boundary  condition,  which, 
n 

for  a  perfect  incompressible  fluid,  is  that  the  velocity  normal  to  the 
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surface,  vn,  oust  vanish.  The  velocity  normal  to  the  surface  Is  given  by 
equation  (3.13)  of  Morrison  et  al4  or 


v 
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rl_  .  (1  -  u2)  dr  L]  , 

L3r  r2  dy  3yJ  f  ' 


(8) 


where  y  *  cos  8  and  r  is  given  by  equation  (6).  Using  equations  (6)  and  (7) 
in  equation  (9)  and  letting  vn  -  0  we  obtain 


...  3Va 
,(1 )  _  _ o 


n 

n+1 


2  5n1 


,  (2) 


An  -’domi  l  *k  C1n,k<1Xln>‘ 


3Va 

8 TTTn  kL  V'  <1k,A>  <£",n>  cu,x(2A  -  <W 


(9) 


where  <ik|n>  are  Clebsch-Gordan  coefficients*7  (In  Rose,7  the 
CUkn;00)  *  <*k|n>),  the  c,  are  given  by 


C  .  -  a(a+l )  ♦  b(b+1 )  -  c(c+1)  , 

ao,c 


and  6  .  is  the  KronecXer  delta  function, 
nl 
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3.  ENERGY 


For  the  change  in  tha  kinetic  energy  of  the  external  fluid,  we  use 


*  •  /  7  PV2  dr  -  /  ^  pV2dr  , 


(10) 


where  the  first  integral  covers  the  region  r(0)  <  r  <  Rq,  and  the  second 
integral  covers  the  region  0  <  r  .<  Rq,  and  in  the  end  result  we  let  Rq  becone 
infinite.  When  equation  (4)  is  used  in  equation  (10),  we  obtain  (RQ  ♦  •) 


T 


(11) 


Thus,  we  need  oily  A1  through  second  order  to  obtain  the  change  in  kinetic 
energy.  All  the  Clebech-Gordan  coefficients  appearing  in  equation  (9)  are 
sieple  for  n  »  1,  and  algebraic  expressions  for  then  can  be  found  in  Rose.7 
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,  (n  -  1)  6n3+  7n2  -  2n  ♦  3) 
n  (2n  -  1 ) ( 2n  +  1)2(2n  +  3) 


_ n  (n  +  1) 

9n  "  (2n  -  1 )  (2n  +  1)(2n  +  3) 


To  obtain  the  result  given  in  equation  (12) ,  we  have  used  the  constraint  of 
constant  volume  as  given  in  Page . 5 , 


a  *  a 


.X  T  an 
"a  *•  2n  +  1 


(13) 


The  surface  energy  of  the  drop  is  s imply  the  product  of  the  surface  tension, 
Y,  and  the  surface  area  of  the  drop  and  is  given  by 


U 

s 


(n  -  1 )  (n  +  2)a 

2*r  I - 
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2n  ♦  1 


(14) 


and  the  constraint  given  in  equation  (13)  was  used  in  Page5  to  obtain  this 
result.  He  now  have  the  total  energy  of  the  system  expressed  in  terns  of  the 


4.  THE  EQUATIONS  FOR  aft 


The  condition  for  mini nun  energy  can  be  written 


J*— 
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(1)  (a  -  ax  )  to  determine  the  shape  of  the  drop  as  a  'function  of  fluid 
n  n 

velocity,  V.  Because  of  the  similarity  of  equation  (18)  to  the  corresponding 
equation  in  Page,5  it  is  expected  that  the  drop  will  deviate  considerably  from 
an  oblate  spheroid  for  moderate  fluid  velocities  or  for  large  drop  radii  a. 
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